Two phase picture in driven polymer translocation 
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Two phase picture is a simple and effective methodology to capture the nonequilibrium dynamics 
of polymer associated with tension propagation. When applying it to the driven translocation 
process, there is a point to be noted, as briefly discussed in our recent article [Phys. Rev. E 85, 
061803 (2012)]. In this article, we address this issue in detail and modify our previous prediction 
[Euro. Phys. J. E 34, 135 (2011)] by adopting an alternative steady-state ansatz. The modified 
scaling prediction turns out to be the same as that of the iso-flux model recently proposed by 
Rowghanian and Grosberg [J. Phys. Chem. B 115, 14127-14135 (2011)]. 
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Tension-propagation mechanism is inherent in the dy- 
namical response of polymers to local stimuli. Sometime 
ago, it has been proposed as an essential ingredient gov- 
erning the driven translocation process [l[. Such a view 
seems to have been well validated in recent studies 0. 

In ref. [3, H| , we analyzed in detail the translocation 
dynamics based on two phase formalism, in which the 
portion of the translocation polymer in cis side is pic- 
tured to be composed of the moving and the yet quies- 
cent equilibrium domains. The point of the formalism is 
to focus on the dynamics of the domain boundary, which 
is a gross variable describing how the driving force at the 
pore site is transmitted along the chain backbone. For 
such a coarse-graining to be meaningful, other degrees 
of freedom need to be properly eliminated, the effect of 
which should be reflected in the behaviors of the gross 
variable. Most naturally, the moving domain is assumed 
to be in a steady-state with a characteristic velocity V(t). 
Then, the steady-state dynamical friction/extension laws 
are coupled to the continuum equation describing the 
conservation of the mass. In our earlier article [3(, we 
set V(t) = Vn(t), where vn(t) represents the segment ve- 
locity around the rear end of the moving domain. This 
steady-state ansatz, however, turned to be not compat- 
ible with the mass conservation relation inherent in the 
translocation process. This point was shortly mentioned 
in the appendix of our article To make the point 
clearer, we present a detailed account in this article, 
where our previous result Q is modified by adopting an 
alternative steady-state ansatz suggested in refs. 
We then point out that the alternative ansatz is inti- 
mately related with the iso-flux model [f|. 

Basic equations — Let us consider the driven transloca- 
tion process as schematically depicted in Fig. 1. A poly- 
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FIG. 1: (Color Online) Illustration of the driven polymer 
translocation process. 



meric chain is initially dissolved in the solution (x < 0), 
taking a random coil conformation with the overall size 
Rq ~ uNq , where Nq is the total segment number, a 
is the segment length, and v is the Flory exponent. As 
those in ref. the length, force and time are made 
dimensionless by the units a, hoTja and r\a? jk^T, re- 
spectively (fceT is the thermal energy; r\ is the viscosity 
of the solvent). 

There is a thin wall at a; = with a small pore, where 
the driving force with the constant magnitude / is locally 
exerted in rr-direction from the cis to the trans side. One 
chain end is sucked at time t = 0. The translocation pro- 
cess, then, proceeds with the tension propagation along 
the chain. Segments are numbered from the first sucked 
end to the other end (A^th segment). According to the 
so-called two phase formalism, the dynamics of driven 
translocation is described as the growth process of the 
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moving tensed domain, which is specified by the segment 
N(t) at the end of the moving domain, the segment M(i) 
at the pore, the size R(t) and the representative (or av- 
erage) velocity V(t) of the moving domain at time t (See 
Fig. 1). These can be determined by the following set of 
equations: 



dN(t) dM(t) 



dt 



dt 



j R {t)-jo{t) 



, . dR . . 

Mt)^, (i) 



V(t)R(t) ~ F" 



N(t) - M(t) ~ Rf-P" (~ Ra ) 



(2) 



(3) 



(4) 



where ^(i) = o- R (t)v R (t) and j (*) = 0"o(*H(*) arc tllG 
segment fluxes at the moving domain boundary and at 
the pore, respectively. As explained in ref. [H, eq. ([1} is 
the continuum equation, eqs. © and ^ are the dynam- 
ical equations of state describing the steady-state rela- 
tion among velocity-extension- force (eq. ([5])) and mass- 
extension-force (eq. ([3])) for a dragged polymer, where 
p z = z — 2 and p v = (1 — v)jv in trumpet regime 
/j</</*, while p z = 1 and p„ = for stronger 
force f*<f 1 - Note that, while the segment line density 
at the pore cro ~ f~ p " is given by the exponent related to 
equilibrium (static) situation, eq. © reflects the dissipa- 
tion mechanism through the exponent p z , which contains 
the dynamical exponent z characterizing the dynamics of 
"critical fluctuation" . 
Ansatz (i) — In ref. Q, we eliminate the pore flux 



dM(t) 
dt 



jo(t) = cr v (t) 



(5) 



in eq. ([T]) and made the steady-state ansatz for the mov- 
ing domain by setting v R (t) = V(t), which leads to 

Combining this with eqs. ([2]), ((4J and the segment line 
density at the boundary a R (t) ~ VR(t)~ q ~ V(t)~ q , 
where q = (1 - i/)/[(z - l)v] for /« < / < /** and g = 
for larger force, one obtains the tension-propagation law 



R(t) * (tfy 



(7) 



with /3 = p z (l — q) and 7 = + v — vq), thus, the 
scaling formula for propagation time is 



K.r 



(8) 



1 For the definition of the characteristic force /j ~ Nq U , /* ~ 1 
and /** ~ IVq, see ref. 0. Strong force regimes are further 
classified as the stem-flower (/* < / < /**) and the strong- 
stretching (/** < /). 



with a = is/"/. 

Ansatz (ii) — An alternative steady-state ansatz is to 
set i>o(i) = V(t) as was done in ref. 0,0] • Then, combin- 
ing eqs. ([3]) and (0, we obtain 



<Jo V(t) 



dR(t) \ 
dt J 



dN(t) 



dt 



(9) 



=> 



fP'-P-Rit)- 1 ~ R{t) {1 - U ^ v 



dR{t) 
dt 



dR(t) 
dt 



(10) 



where eqs. ([2]), (|4]) and the expression for the line density 
at the pore (see eq. ©) are used. The tension propaga- 
tion law and the propagation time in leading order can 
be expressed in the form of eqs. (O and ©, respectively, 
but now the exponent (3 and 7 are modified as j3 = p z —p v 
and 7 = u/(l + v). The deduced exponents coincide with 
those based on the iso-flux model [all] (see below). 

The difference in the dynamical evolutions of the num- 
ber M(t) of the translocated segments depending on the 
adopted ansatz is plotted in Fig. [5] The inflection points 
in Fig. [2ja) correspond to the end of the propagation 
stage, when most of segments have been already sucked. 
The subsequent post-propagation stage can be analyzed 
using eqs. ©/© 0, where, in contrast to the propaga- 
tion stage, the growth rate in M(t) increases with time 
due to the decrease in the overall friction with the process 
advanced. As shown in Fig. 2(b), the growth of M(t) in 
the propagation stage follows the asymptotic dynamical 
scaling M(t) ~ t 7 /^, with the slight deviation in the early 
period due to the finite-size effect. 

Iso-flux model- We now attempt to resolve the dis- 
crepancy between the above two results. The analysis 
of eqs. ([1]) , and (@| reveals a unique property of the 
moving domain in the translocation process. Given the 
average relation eq. (0| between R(t) and N(t) , one can 
show that the last terms in eqs. |T]) and ((3|) are irrele- 
vant, i.e., dN(t)/dt > cr R (t)dR/dt and N(t) > Rf~ p » 2 . 
It then follows that 



and 



M{t) ~ N(t) 



j R (t)~j (t) 



111) 



(12) 



In the scaling level, the fluxes into and out of moving 
domains balance, as was recently proposed by Rowgha- 
nian and Grosberg To check the validity of the 



First inequality: Let us write the force balance on the rear end 
blob v R ~ Since £ R < R, this leads to the inequal- 

ity R~~ 2 v R > 1/R. Eliminating v R through a R ~ g R /£ R ~ 
v^ 1 i)"^ onc g n( j s j\r 1— " > cr R , which is equivalent to 

be proven with eq. Q- Second inequality: Using eq. Q, this 
can be rewritten as fO-~"^l u rO-—")/" > 1, which holds under 
the condition fR > 1 of interest. 
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FIG. 2: (Color Online) Number M(t) of the translocated 
segments as a function of time calculated according to the 
steady-state ansatz (i) or (ii) with v = 3/5, z = 3 and No — 
100. (a) Entire process in normal scale, (b) Propagation 
stage in double-logarithmic scale. The dashed lines represent 
the asymptotic scaling with the slope ~f/v. 



steady-state approximation, let us set V(t) = v(t,x) 
with 3x € [-R(t),0] so that v R (t) = V(t) - 5v R {t) and 
v (t) = V(t) + Sv (t). The iso-fiux condition eq. $TZ§ is 
written as 



a [V(t) + Sv (t)] 



<TR[V(t) 

Svr 
V 



1 - 



- Sv R (t)} 

1-8 



(13) 



where a R (t)/a ~ (1 - Sv R /V)- q W(t) q with W{t) = 
fR(t) for f s < f < f* and W{t) = f*R*(t) for /* < 
/ 3 . Given the longest relaxation time r rc i ax ~ R/V 
of the steadily driven domain and the typical shear rate 
7 ~ 5v/R = (vq—v r )/R, the steady-state approximation 
implicitly assumes the condition r rc i ax 7 < 1 Sv(t) < 
V(t). This self-consistency condition is satisfied for the 
ansatz (ii), i.e., 8v a = 6v R = [1 - W(*) _9/(1 ~ ?) ]V < 
V, but not for the ansatz (i), i.e., Sv R = => Svq = 
[W(t) q - 1}V in the force range /((</< /** where 
W{t) q > 1. 



To summarize, the unique driving mode inherent in the 
translocation process is manifested in the iso-fiux con- 
dition, which is to be contrasted with the mechanical 
driving by pulling chain end 0, 0, §■ In the latter, the 
moving domain is constituted of all the N(t) segments, so 
that eqs. ([3]) and ((H) are modified as N(t) ~ Rf~ Pv and 
Nty)" <C R(t), respectively. It has turned out that our 
previous ansatz and corresponding predictions Q are in- 
compatible with the iso-fiux condition. From purely the- 
oretical point of view 0-Q , the scaling of the propagation 
time (thus, the translocation time) is given by 



n] +1 



f- 



(1< 



f) 



/<!) 



(14) 



It is somewhat puzzling that the prediction based on the 
ansatz (i) seem to be better correlated with numerical 
and real experiments 0| ■ This might be ascribed to finite- 
size effects and various crossovers detailed in ref. 0]; (i) 
The last term in eq. (|10[) serves as a correction term to 
the tension-propagation scaling law (see the inequality 
above eq. (jlip ). (ii) The post-propagation stage adds a 
correction term to the translocation time, (iii) Iso-fiux 
condition (cqs. (TIT]) and (fT12j) ) is also influenced by the 
finite-size correction. These effects may be enhanced in 
realistic situations due to the wall and/or the trans-side 
effects, which would be dependent on the size and proper- 
ties of the pore etc. Since the pore flux would be sensitive 
to such specific pore effects, one may then be tempted to 
remove the pore flux (eq. (O) from eq. (fT]) and to follow 
the route of the ansatz (i). If the outgoing flux jo(t) is 
suppressed enough compared to that expected from the 
iso-fiux condition, the description based on the steady- 
state ansatz (i) may attain its apparent suitability. A 
very recent theoretical study based on a more "accurate" 
two phase description seems to support such a view 0], 
where it has been shown that the frictional interaction 
between the pore and the polymer significantly affects 
the dynamics for polymers with experimentally relevant, 
finite chain length. 
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